Evolution of a Schwarzschild Black Hole in Phantom-like 

Chaplygin gas Cosmologies 



Mubasher Jamir 

Center for Advanced Mathematics and Physics, 
National University of Sciences and Technology, 
E&ME campus, Peshawar road, Rawalpindi - 46000, Pakistan 

July 31, 2009 



Abstract 

In the classical relativistic regime, the accretion of phantom energy onto a black hole 
reduces the mass of the black hole. In this context, we have investigated the evolution of a 
Schwarzschild black hole in the standard model of cosmology using the phantom-like modi- 
fied variable Chaplygin gas and the viscous generalized Chaplygin gas. The corresponding 
expressions for accretion time scale and evolution of mass have been derived. Our results 
indicate that mass of the black hole will decrease if the accreting phantom Chaplygin gas vi- 
olates the dominant energy condition and will increase in the opposite case. Thus our results 
are in agreement with the results of Babichev et al 21j who first proposed this scenario. 
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1 Introduction 



The evidence of accelerated expansion in the observable universe is quite compelling and has been 
confirmed by various astrophysical investigations including observations of supernovae of type la 
[H [2] , anisotropics of the cosmic microwave background radiation [3l 0] , large scale structure and 
galaxy distribution surveys [5]. This expansion of the universe is supposedly driven by an exotic 
energy commonly called 'dark energy' possessing negative pressure p < and positive energy 
density p > 0, related by equation of state (EoS) p = top. It should be noted that p = up is not a 
true EoS for dark energy rather a phenomenological description valid for a certain configuration 
[6]. Astrophysical data suggests that about two third of the critical energy density is stored in 
the dark energy component. The corresponding parameter lo is then constrained in the range 
— 1.38 < uj < —0.82 [7j. It shows that the EoS of cosmic fluids is not exactly determined. The 
genesis of this exotic energy is still unknown. The simplest and the earliest explanation of this 
phenomenon was provided by the general theory of relativity through the cosmological constant 
A. The observational value of its energy density is 56 to 120 orders of magnitude smaller than 
that derived from the standard theory [8]. The satisfactory explanation of this phenomenon 
requires extreme fine tuning of the cosmological parameters. Other problem associated with A 
is the coincidence problem (i.e. Why did the cosmic accelerated expansion start in the presence 
of intelligent beings? alternatively why the energy densities of matter and dark energy are of the 
same order at current time?) which is as yet explained either through anthropic principle [9], 
variable cosmological constant scenario |10j or by invoking a dark matter-dark energy interaction 
[m [121 [13] . Ill this context, several other models have been proposed among them are models 
based on homogeneous and time dependent scalar fields termed as quintessence [H], quintom 
[151 [16] aiid k-essence [T7], to name a few. 

The interest in phantom energy arose when Caldwell et al [18] explored the cosmological 
consequences of the EoS, uj < —1. The dark energy can achieve this EoS if it is assumed to be 
variable quantity i.e. u;(z), where z is the redshift parameter. Thus u) evolves as: for matter 
dominated universe a; = 0, in quintessence phase —1 < w < —1/3, for cosmological constant 
dominated arena a; = — 1, while in phantom regime a; < — 1. This scenario appears to be 
consistent with the observations [19]. In phantom cosmology, the energy density of phantom 
energy will become infinite in a finite time leading the 'big rip', a kind of future singularity. 
Moreover, due to strong negative pressure of phantom energy, all stable gravitationally bound 
objects will be dissociated near the big rip. These findings were later confirmed in [20] by doing 
numerical analysis for the solar system and the Milky Way galaxy. In this context, the accretion 
of phantom dark energy onto a black hole was first modeled by Babichev et al ^21j who proved 
that black hole mass will gradually decrease due to strong negative pressure of phantom energy 
and tend to zero near the big rip where it will finally disappear. Note that lo > —1 leads to the 
opposite scenario where black hole mass increases by accreting dark energy until its event horizon 
swells up to swallow the whole universe |22j . Later studies [23] showed that quantum effects 
dominate near the big rip singularity and consequently the mass of black hole although decreases 
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but stops decreasing at a finite value. In another investigation [24j, it was demonstrated that 
the physical black hole mass will increase due to accretion of phantom energy consequently the 
black hole horizon and the cosmological horizon will coincide leading the black hole singularity to 
become naked, all in a finite time. This analysis has been extended for the Riessner-Nordstrom, 
Kerr-Neumann and Schwarzschild-de Sitter black holes as well [251 ESJ EH ESI ESI [30] . This result 
apparently refutes the cosmic censorship conjecture (or hypothesis) which forbids the occurrence 
of any naked singularity. However, the formation of naked singularities is not completely ruled 
out. Numerical simulations of gravitational collapse of spheroids show that if the collapsing 
spheroid is sufficiently compact, the singularities are hidden inside black hole while they become 
naked (devoid of apparent horizon) if the spheroid is sufficiently large [31j . The future singularity 
of big rip is alternatively avoided by the 'big trip' where the accretion of phantom energy onto 
a wormhole will increase the size of its throat so much to engulf the whole universe [32[ [33] . 
Another interesting scenario appears in cyclic cosmology where the masses of black holes first 
decrease and then increase through the phantom energy accretion and never vanish [34] . The 
implications of generalized second law of thermodynamics to the phantom energy accretion onto 
a black hole imply that accretion will be significant only near the big rip. If this law is violated 
than the black hole mass will decrease |35] . The thermodynamical investigations of phantom 
energy imply its positive definite entropy which tends to become constant if the phantom energy 
largely dominates the universe [36]. This results in the late universe to be hotter compared to 
the present. 

We here discuss the accretion of phantom like modified variable Chaplygin gas and the viscous 
Chaplygin gas seperatly onto a black hole. This accretion of the phantom fluid reduces the mass 
of the black hole. This works serves as the generalization of the earlier work by Babichev et al 
\21\ [37] who initiated the concept of accretion of exotic matter on the black hole. We have built 
our model on the same pattern by choosing more general EoS for the dark energy. 

The outline of the paper is as follows: In the next section, we discuss the relativistic model of 
accretion onto a black hole. In third section, we investigate the evolution of the mass of black 
hole by the accretion of modified variable Chaplygin Gas (MVG) while in the fourth section, 
we discuss the similar scenario with the viscous generalized Chaplygin gas (VCG). Finally, we 
present conclusion of our paper. The formalism adopted here is from [21] . 

2 Accretion onto black hole 

We consider a Schwarzschild black hole of mass M which is gravitationally isolated and is 
specified by the line element (in geometrical units c = 1 = G) : 




(1) 
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The black hole is accreting Chaplygin gas which is assumed to be a perfect fluid specified by 
the stress energy tensor 

T^''' = {p + p)ui'u'' -pg^'''. (2) 

Here p and p are the pressure and energy density of the Chaplygin gas respectively. Due to 
static and spherically symmetric nature of the black hole we assume the velocity four vector 
= (u*(r), 77,''(r), 0, 0) which satisfies the normalization condition it'^it^ = —1. Thus we are 
considering only radial in-fall of the Chaplygin gas on the event horizon. Using the energy- 
momentum conservation for T'*'^, we get 



ux\p + p)\ l --+V? = Ci, 



X 



(3) 



where x = r/M and u = vJ" = dr / ds \s the radial component of the velocity four vector and 
Ci is a constant of integration. The second constant of motion is obtained from u^Tf^^ = 0, 
which gives 

r Ph 

dp' 



ux exp 



pf+p'ip') 



(4) 



where A \s a, constant of integration which is determined below for two models of Chaplygin 
gas. The quantities poo and ph are the densities of Chaplygin gas at infinity and at the black 
hole horizon respectively. Further using Eqs. (3) and (4), we obtain 



I '''' / 

(p + p)y'l-^ + n2exp -jj^ 



C2, 



(5) 



where C2 = — Ci/A = p^o +p{poo)- In order to calculate M, the rate of change of mass of 
black hole we integrate the Chaplygin gas flux over the entire horizon as, M = / T^dS where 
Tl denotes the momentum density of Chaplygin gas in the radial direction and dS = ^/^dOd^p 
is the surface element of black hole horizon. Using Eqs. (2-5), we get this rate of change as 



— = AitAM^{p + p). 



Integration of Eq. (6) yields 



M = M,- 



T 



(6) 



(7) 



which determines the evolution of mass of black hole of initial mass Mj and r is the characteristic 
accretion time scale given by 

(8) 



r ^ = 4TrAMi{p + p). 
The number density and energy density of Chaplygin gas are related as 



njph) 
n{poo) 



exp 



Ph 



dp' 

7+7(7) 



U'oo 



(9) 
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here n{ph) and n(/>oo) are the number densities of the Chaplygin gas at the horizon and at 
infinity respectively. Further the constant A appearing in Eq. (8) is determined as 



n{ph) 2 



ux^ = -A, (10) 



which is an alternative form of energy momentum conservation Eq. (4). Moreover, the critical 
points of accretion (i.e. the points where the speed of flow achieves the speed of sound V"^ = 
= dp /dp) are determined as follows 
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= — V"^ = fill 

where V"^ = ~ 1- Finally, the above Eqs. (9 - 11) are combined in a single expression 

as 

P*+P*(P*) _[,^3^2(^^)]l/2^^oo+p(Poo)_ (^2) 



n(p*) n{pc 



3 Accretion of modified variable Chaplygin gas 



The Chaplygin gas had been proposed to explain the accelerated expansion of the universe [38] . 
It is represented by a simple EoS, p = —A'/p where A' is positive constant. The corresponding 
expression for the evolution of energy density is 

P = y^' + l' (13) 

where S is a constant of integration. From Eq. (13) we find the following asymptotic behavior 
for the density |39j : 

p ~ VBa-^, a « (B/A')^/^ (14) 

pr^pr^VA', a^{B/A'f'^. (15) 

Thus for small a, it gives matter dominated era at earlier times while for large a we get dark 
energy dominated era at late times. Thus Chaplygin gas has the property of giving a unified pic- 
ture of the evolution of the universe. The observational evidence in favor of cosmological models 
based on Chaplygin gas is quite encouraging [iOl HH [42l [43] . The Chaplygin gas model favors 
a spatially fiat universe which agrees with the observational data of Sloan Digital Sky Survey 
(SDSS) and Supernova Legacy Survey (SNLS) with 95.4 % confidence level [44]. Consequently, 
various generalizations of Chaplygin gas have been proposed in the literature to incorporate any 
other dark component in the universe (see e.g. |45 t 1461 Wh I48j and references therein). 

We here consider an equation of state which combines various EoS of Chaplygin gas given by 

m 

p = A'p-^, B{a) = B,a-^. (16) 
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Here A', Bq and ra arc constant parameters with < a < 1. For A' = 0, Eq. (16) gives 
generalized Chaplygin gas. Further \i B = and a = 1, it yields the usual Chaplygin gas. 
Also Eq. (16) reduces to modified Chaplygin gas if only B = Bo- Moreover, if only A' = 0, the 
same equation represents variable Chaplygin gas. 

We consider the background spacetime to be spatially fiat {k = 0), homogeneous and isotropic 
represented by Friedmann-Robertson- Walker (FRW) metric. The spacetime is assumed to con- 
tain only one component fluid i.e. the phantom energy represented by the Chaplygin gas EoS. 
The corresponding field equations are 



(17) 



where = Stt/S. The conservation of energy is 

p + 3i7(p + p) =0. 
Using Eqs. (16) and (19), the density evolution is given by 

3I?„(1 + a) 



P 



+ 



+ a){l + A') - m} a3a+")(i+^')_ 



(18) 



(19) 



(20) 



Here ^' is a constant of integration. Note that to obtain the increasing energy density of phantom 
energy with respect to scale factor a{t), we require the coefficients of a{t) in Eq. (20) to be 
positive i.e. ^ > 0, Bo{l + a) > and 3(1 + a)(l + A') — m > 0. Moreover, the exponents 
of a{t) must be negative i.e. m < and 3(1 + a)(l + A') < to obtain increasing p. These 
constraints together imply that m > 3(1 + a)(l + A'). Another way of getting positive p is by 
setting m > 0, 1 + ^' > and m < 3(1 + a)(l + A'). Further, using Eq. (9) the ratio of the 
number density of Chaplygin gas near horizon and at infinity is calculated to be 



njph) 
n(poo) 



A') - B(a) 



p^"(l + A') 



B(a) 



(l+a)(l+A') 



(21) 



Notice that the function B(a) can be expressed in terms of p implicitly and is determined from 
Eq. (20). Making use of Eq. (11), the critical points of accretion are given by 



I + 3A2 



I + 3A2 
2A2 



where 



., aB(a) . 



(22) 



(23) 



Thus for the critical points to be finite and positive, we require either A2 > or A2 < and A2 < 
— 1/3. For the accretion to be critical, the quantity must become supersonic from the initial 
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subsonic somewhere near the black hole horizon. For the MVG, we have uj = A' - B / p^+'^ < 0, 
since A' < —1. One can observe that fluids having EoS w < are hydrodynamically unstable i.e. 
the speed of sound in that medium can not be defined since Cg < 0. In order to overcome this 
problem Babichev et al [50j proposed a redefinition of w with the help of a generalized linear 
EoS given by p = (3{p — po), where f3 and po are constant parameters. Here /? > refers to a 
hydrodynamically stable while /? < corresponds to hydrodynamically unstable fluid. We will 
not be interested in the later case here. Note that now two parameters uj and /3 are related by 
u) = P{p — Po)/ P- Further lij < now corresponds to /? > and po > p thereby making the 
previously unstable fluid, now stable. We also have = dp/ dp = (3. Since for stability, we 
require /? > and < < 1, it leads to < j^{A' p - B / p'^) < 1 and < /3 < 1. Hence the 
EoS parameter is now well-defined with A' < —1 and po > p. Thus the stability of the phantom 
like MVG is guaranteed with the use of generalized linear EoS. 

The constant A is determined from Eq. (10) to give 



4 V A2 , 

Using Eq. (8) the characteristic evolution time scale becomes 



.-=.A,,p + p)^(i^)'^ (25) 



Using Eq. (25) in (7), the black hole mass is given by 



M{t) = Mi 



l-T:Mit{p + p)- ^ 



Al /1 + 3A2\ 



3/2' 



As / 



(26) 



which determines the evolution of mass of black hole accreting phantom MVG. It can be seen 
that if the phantom MVG violates the dominant energy condition p + p > than mass M of 
the black hole will decrease. Contrary if this condition is satisfied than M will increase. Thus 
in the classical relativistic regime, this result is in conformity with the result of Babichev et al 
|21j . We also stress here that although our metric (1) is static, we get a dynamical mass M{t) 
in Eq. (26). Astrophysically the mass of a black hole is a dynamical quantity: the mass will 
increase if the black hole accretes classical matter (which satisfies p + p > 0) however it will 
decrease for the exotic phantom energy accretion. The mass can also decrease if the Hawking 
evaporation process is invoked. Hence the static black holes may not necessarily correspond to 
the astrophysical black holes. We also stress that cj > (cj < 0) corresponds to non-phantom 
(phantom) MVG fluid; although the accretion through the critical point is possible in both the 
cases, only phantom MVG violating the dominant energy condition will reduce the mass of black 
hole. 
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4 Accretion of viscous generalized Chaplygin gas 



In viscous cosmology, the presence of viscosity corresponds to space isotropy and hence is im- 
portant in the background of FRW spacetime [51^ [52l I53j . The presence of viscous fluid can 
explain the observed high entropy per baryon ratio in the universe |54) . It can cause exponential 
expansion of the universe and can rule out the initial singularity which mares the standard big 
bang picture. The matter power spectrum in bulk viscous cosmology is also well behaved as 
there are no instabilities or oscillations on small perturbation scale j55j . Any cosmic fluid having 
non-zero bulk viscosities has the EoS Pe// = P + n, where p is the usual isotropic pressure and 
n is the bulk viscous stress given by 11 = £,{p)u^^ [56]. The scaling of viscosity coefficient is 
C = Co/o" where n is a constant parameter and i{to) = Co- Note that for < n < 1/2, we 
have de Sitter solution and for n > 1/2 we get deflationary solutions. The viscosity coefficient 
is generally taken to be positive for positive entropy production in conformity with the second 
law of thermodynamics [57]. Moreover, the entropy corresponding to viscous cosmology is al- 
ways positive and increasing which is consistent with the thermodynamic arrow of time. Infact 
the cosmological model with viscosity is consistent with the observational SN la data at lower 
redshifts while it mimics the ACDM model in the later cosmic evolution [58]. It is proved in 
|59t [60] that FRW spacetime filled with perfect fluid and the bulk viscous stresses will violate 
the dominant energy condition. 

Thus the effective pressure is given by 



where 11 = —3H^ and p = x/p°' with x is a constant. Thus in the VCG case, the standard FRW 
equation becomes [HI] 



which shows that the viscosity term serves as the source term. Using Eqs. (17) and (27) in (29), 
we get 



Peff = p + n. 



(27) 




(28) 



Further the energy conservation principle gives 



p + 3H{p+Peff)=0, 



(29) 




3<(p)^ = 0. 



(30) 



Thus solving Eq. (30) we get 



p 



3 




(31) 



For our further analysis we shall assume ^ to be constant. 
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The ratio of the number density of VCG near black hole horizon and at infinity is given by 



njph) 



exp 



Ph 

/■ 



dp' 



= A, 



(32) 



The corresponding critical points of accretion are 

A4 



3A4 



1 



3A4 - 1 
2A4 '' 



where 



= A4 



(33) 



(34) 



.P* ' * 

Notice that for the critical points to be finite and positive valued we require either A4 < or 
A4 > 1/3. Using Eq. (11) we see that the speed of flow at the critical point is V'^ = — A4. 
Further, the EoS parameter is u; = x/P^"*"" ~ ^^H/p (= x/P^"*"" ~ V^'t^/y^)- Note that if 
X < then a; < and stability of VCG is lost. However, if we here invoke the argument 
presented in the last section, we can consider accretion with a; < 0. Using the generalized 
linear EoS p = (3{p — po) for the phantom energy, we obtain /3 > and Po > p for u < 0. 
Using the definition = dp/dp = P and stability requirements /? > and < < 1 lead to 
< pj^p ix/p" ~ \/3/JK^) < 1 and < /3 < 1. The EoS parameter (3 is now well-defined with 
X < and po > p- Therefore the stability of the phantom like VCG is assured with the use of 
generalized linear EoS. 



Using Eq. (10) the constant A is now determined to be 

3A4-n^/^ 



-A = A3 

The characteristic evolution time scale is 



2A4 



r"^ = 47rMi(p + p)A3 



3A4 - 1 
2A4 



3/2 



(35) 



(36) 



Using Eqs. (35) and (36) in (7), we get black hole mass evolution as 



M{t) = Mi 



1 - 4TTMit{p + p) A3 



3A4 - 1 
2A4 



3/2' 



(37) 



It can be seen that black hole mass will decrease when p + p < and increase in the opposite 
case. It is emphasized that this result is valid till the contribution of viscous stress is negligible 
compared to isotropic stress. For the sake of clarity, we emphasis that fluid violating the standard 
energy conditions is termed 'exotic' and hydrodynamically unstable i.e. its existence is not fully 
guaranteed. But this conclusion is drawn due to the 'bad' choice of the EoS {p = up) in the 
analysis. The result is reversed and remedied when we introduce the generalized linear EoS in 
our model which makes the accretion of exotic fluid much more practical. 
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5 Conclusion 



We have investigated the accretion of two different forms of phantom-hke Chaplygin gas onto a 
Schwarzschild black hole. The time scale of accretion and the evolution of mass of black hole are 
derived in the context of two widely studied Chaplygin gas models namely the modified variable 
Chaplygin gas and the viscous generalized Chaplygin gas. Although the phantom energy is an 
unstable fluid as it corresponds to a medium with indeterminate speed of sound and super- 
luminal speeds. These pathologies arise due to bad choices of the equations of state for the 
phantom energy and hence can be removed by choosing some suitable transformation from one 
EoS to another or a totally new EoS for this purpose. This work serves as the generalization of the 
earlier work by Babichev et al [21]. It should be noted that we have ignored matter component 
in the accretion model. Thus it will be more insightful to incorporate the contributions of matter 
along with the Chaplygin gas during accretion onto black hole. Moreover our analysis can be 
extended to the case of rotating black holes as well. 
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